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' Abstract. It is shown how the strong interaction dynamics of a multi-quark Lagrangian 

CS| , affects the catalysis of dynamical symmetry breaking by a constant magnetic field in (3 + 1) 

dimensions. Attention is drawn to the local minima structure of the theory. 
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^ '. 1 Introduction 

m 

0^ ■ The existence of a zero-energy surface in the spectrum of a Dirac particle is ensured for any 

homogeneous magnetic field with a fixed direction by a quantum mechanical supersymmetry of 
the corresponding second-order Dirac Hamiltonian |T]. For interacting fermions subject to an 
• attractive two-body force in 2 -|- 1 and 3-1-1 dimensions, a constant magnetic field catalyzes 

OO . the dynamical symmetry breaking leading to fermion pairing and to a fermion mass even for 

arbitrarily weak coupling strength [21 [31 14] . The zero-energy surface of the lowest Landau level 
(LLL) plays a crucial role in the dynamics of such fermion pairing [3 [H [3, and the generated 
^ . fermion mass, M^yn, turns out to be much smaller than the Landau gap ~ y^\eH\. The dynamics 

^ I of the fermion pairing in the homogeneous magnetic field is essentially (1 -|- l)-dimensional, 

and the deep analogy of this phenomenon with the dynamics of electron pairing in BCS [8] 
has been stressed. The generation of mass via catalysis does however not prevail for arbitrary 
magnetic field configurations. For instance, it has been demonstrated in [9l[T0] that the Nambu- 
Jona-Lasinio (NJL) model [11] minimally coupled to a background magnetic field with variable 
direction becomes massive only beyond some critical value of the strong coupling constant. In 
these cases dynamical chiral symmetry breaking corresponds to the standard scenario. 

In the present paper we show and discuss the influence that the structure of interaction 
terms, specified below, may have on the effect of catalysis by a constant magnetic field (from 
now on we use interaction terms to denote interaction among the fermions, the interaction with 
the magnetic field will be explicitly referred as such). We obtain that catalysis of dynamical 
symmetry breaking always occurs at arbitrarily small magnetic fields for a system which is 
subcritical in the interaction couplings in absence of the magnetic field. However for increasing 
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strength of the magnetic field and depending on the values of interaction couplings, another 
minimum of the effective potential emerges at larger dynamical fermion masses, which competes 
with the catalyzed minimum and eventually becomes the globally stable state of the system at 
characteristic scales H = 10^^ Gauss. The broken symmetry is not restored with increasing 
strength of the magnetic field. 

A few words about the fermion interaction model used in this paper: in recent years we 
have reported on the importance of adding eight-quark interaction terms |12] to the four-quark 
f^(3)L X U{2>)r chiral symmetric NJL Lagrangian joined with the U{1)a breaking 't Hooft six- 
quark interaction [13]. Their inclusion leads to a scalar effective potential which is globally 
stable and to a full equivalence between the stationary phase approach (SPA) and the mean 
field (MF) approximation methods used to obtain it. Without them the model is afflicted by 
an unstable (SPA) respectively metastable (MF) vacuum [H], as presented two years ago at 
the Kyiv 2005 conference [15]. There we have also shown that going beyond the leading order 
results of SPA did not cure the instability of the effective potential. Although we prefer to 
consider the inclusion of eight quark interactions as the minimal chiral approach to render the 
vacuum of the NJL model with 't Hooft interaction stable, QCD regarded as an effective theory 
displays an infinite chain of multiquark interactions [16j . The issue of the corresponding Nc 
counting rules is however still open. In the framework of the instanton gas model beyond the 
zero mode approximation multiquark interactions are delivered at any 2n-quark order [n > 2) 
with equal weight in the large N^. limit [1^. On the other hand lattice calculations of gluon 
correlators [18] give a larger weight to the lower correlators. This leads to the expectation 
that after integrating out the gluonic degrees of freedom, a similar hierarchy could arise for the 
multiquark interactions. In our approach higher order interactions are assumed to be suppressed 
in the large Nc counting [121 IE] . 

Since then the eight-quark extended model has been carefully studied in connection with the 
low lying spectra and characteristics of the pseudoscalar and scalar nonets [19], temperature 
dependence of chiral transitions [20], and [21] for the two flavor case, and coupling to a constant 
magnetic field [22]. We have learned that the effects of the eight quark interactions on the mass 
spectra are not essential, except for the mass of the lowest state of the mixed singlet-octet scalar 
states, which decreases with increasing strength of the eight quark interaction terms that vio- 
late the Okubo-Zweig-Iizuka (OZI) rule [23]. However, their effect can be dramatic on studies 
involving changes of the extrema of the effective potential, which are for instance the cases re- 
lated with the thermodynamics, the in medium dependence, and the coupling to electromagnetic 
fields of the multi-fermion system. The latter is the main subject of the present contribution. 
This work is very much based on our recent paper [22], it is however a more extensive version. 
It contains also truly novel calculations, as we chose as starting configurations a regime of com- 
pletely sub-critical solutions (i.e. only the trivial vacuum with vanishing condensate exists in 
absence of the magnetic field), which were not analyzed before for the enlarged multi-quark 
Lagrangian. We study then how the system responds to coupling to a magnetic field, starting 
with an arbitrarily weak strength and increasing it successively. First appear the well known 
catalyzed solutions as minima of the gap equation. Then the fate of the system may be either 
a first order transition or a crossover to heavier condensates, depending on the strength of the 
couplings of multiquark interactions. The paper is structured as follows. In Section[2]we present 
the model Lagrangian in presence of a constant magnetic field and the corresponding effective 
potential. In Section [3] we obtain the solutions of the gap equations and their relation with the 
Landau levels and highlight the differences to the well-known cases. In Section H] we present the 
conclusions. 
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2 Model Lagrangian 

In the hadronic QCD vacuum the relevant fermionic interactions can be modeled by the well- 
known Nambu-Jona-Lasinio (NJL) Lagrangian [11] . The coupling to a constant magnetic field 
has been done in [6] for one-flavour fermions. 

For a more realistic analysis one should take into account the different quark flavours of the 
QCD vacuum. These extensions of the NJL model are well-known [241 125^ [26]. for instance, the 
four-quark U{3)lx U{3)r chiral symmetric Lagrangian together with the U{1)a_ breaking 't Hooft 
six-quark interactions has been extensively studied at the mean-field level [271 [Ml EH [30l [31] . As 
mentioned in the introduction, it has been recently shown [12] that the eight-quark interactions 
are of vital importance to stabilize the multi-quark vacuum. 

The multi-quark dynamics of the extended NJL model is described by the Lagrangian density 

>Ccff = q{i^^D^ - m)q + dg + C^q + Csq, 

where the gauge covariant derivative is equal to = + iQA^ with the external elec- 
tromagnetic field Af^ and quark charges Q = e ■ diag(2/3, —1/3, —1/3). We consider a constant 
magnetic field, = —Hy, Ay = A^ = (Landau gauge). It is assumed that quark fields have 
colour (Nc = 3) and flavour (iVj = 3) indices. The current quark mass, m, is a diagonal ma- 
trix with elements diag(mu, m^, rrig), which explicitly breaks the global chiral ^[/^(S) x SUr{3) 
symmetry of the Lagrangian. We shall neglect this effect in the following assuming that m = 0. 

The most general form of multi-quark interactions in the scalar and pseudoscalar channels 
are of the form 

= [{qKqf + [qijbKqf] , Cqq = K{det qPiq + det qPnq), Csq = Cg^ + Cf^. 

The C/(3) flavour matrices Aq, a = 0, 1, ... ,8, are normalized such that tr(AaAb) = 25ab- The 
matrices Pl,r = (1 =F 75)/2 are chiral projectors and the determinant is over flavour indices, 
which are suppressed here. The determinantal interaction breaks explicitly the axial U{\)a 
symmetry |13j and Zweig's rule. The eight-quark spin zero interactions are given by 

(1) 2 (2) 

= ^91 [{qiPRqm){qmPLqi)\ , C\q = i^g2{qiPRqm){qmPLqj){qj PRqk){qkPLqi) ■ 

G, K, gi, g2 are dimensionful coupling constants: [G] = M~^, [k] = M~^, [gi] = [^2] = in 
units h = c = 1. The term proportional to gi violates the OZI rule. 



2.1 Effective potential 

We proceed by calculating the effective potential of the theory, y(m„, m^, rris), in the constant 
magnetic field Ax = —Hy, Ay = A^ = 0. The arguments, mj, are simply real parameters; they 
are not to be identified with the masses of any presumed one-particle states. Instead, we shall 
use the capital letter Mj for the point where V takes its local minimum, which specifies the 
masses of constituent quark fields. 

The potential is built of the following two terms 

Viniu^ind^ms) = Vst + Vs- (1) 

The first contribution results from the many-fermion vertices of Lagrangian £cff ) after reducing 
them to a bilinear form with help of bosonic auxiliary fields, and subsequent integration over 
these fields, using the stationary phase method. The specific details of these calculations and 
the result are given in [1^. We obtain 

V.t = ^ Uchf + Khuh^hs + ^ {h^f + 3g2hj] , 
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where /i? = + h'^ + h"^, and hf = hf^ + h"^ + . The functions hi depend on the coupHng 
constants G, k, gi, 52 and on the independent variables Aj = rrii — rhi. To find this dependence 
one should solve the system of cubic equations 



(2) 



Ghu + 


A„ + 


lb 


— Khi + 




= 


Ghd + 


Ad + 


—Khs + 
lb 


-jhdh^ + 




= 0, 


Ghs + 


A, + 


T^Khd + 
lb 


-^h,hi + 


2 ' 


= 0. 



In the parameter range 

91 > 0, gi + 3g2 > 0, 



91 V16/ 



(3) 



the system has only one set of real solutions, and this guarantees the vacuum state of the theory 
to be stable [12]. 

The second term on the r.h.s. of equation ([1]) derives from the integration over the quark 
bilinears in the functional integral of the theory in presence of a constant magnetic field H. As 
has been calculated by Schwinger a long time ago |32j 



ys= Yl ysim^,\Q^H\), 



where 



00 

Vs{m,\QH\) = ^ I ^e-'^'p{s,A^)\QH\cotHs\QH\), 



up to an unessential constant. Here the cutoff A has been introduced by subtracting off suitable 
counterterms to regularize the integral at the lower limit, i.e., p{s, A^) = 1 — (1 + sA^)e~*^^ . For 
the fermion tadpole this works as the four-momentum covariant cutoff in the Euclidean space: 
+ P4 < A^. As a result we obtain 



Vsim,\QH\) 



87r2 



A^\QH\ 



+ ^(ln- 



d_ 

dv 



In 27r - 2 In r 



A^ + m? 



2\QH\ 



+ m'^\QH\ In 1 + 



A2 



h? -t- m? 



A2 



2\QH\ 



u=0 



2\QH\ 

The quantity Ci'^^x) denotes the generalized Riemann zeta function [33j . 



3 Analysis of solutions 

In this section we study the extrema of the effective potential, given by the gap equations. We 
start by considering the known case with four-quark interactions in Subsection 13.11 A detailed 
study of several relevant limits and the role of Landau levels in the catalysis of dynamical 
symmetry breaking is presented. Then in Subsection 13.21 the extended system with six and eight 
quark interactions is analyzed and its new features discussed. 
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3.1 The standard four-quark system in a constant magnetic field 



We consider first the simple SU (3) flavour limit for the situation in which m = and n = 
Qi = 92 = 0. For purposes of illustration, we ignore in the remaining the charge difference of u 
and d, s quarks. The averaged common charge |Q| = |4e/9| will be used. In this case one gets 
the potential V{ra) = Nf{m'^/4G + Vs{m, \QH\))}] and the gap equation is obtained as 



dVjm) 
dm 



1 



mNcNfs 

^ 1 2GNr 47r2 



+ 



\QH\ 
47r2 



2\QH\ 



2\QH\ 



^ ' 2\QH] 



0, 



where 'ip{x) = dlnT{x)/dx is the Euler dilogarithmic function, and 
Ji(m2) = ln( 1 + ^ 



A2\ A2 



A^ + 



The gap equation always has a trivial solution, m 
in the equation 



0. The nontrivial solution is contained 



f{m^;A,\QH\)^i^ 



A^ + 
2\QH\ 



\QH\ 
A2 



Ji{m^)+2ln- 



A2- 



2\QH\ 



2\QH\ 



(4) 



Note that although our regularization scheme differs essentially from the one used in [6], the 
main conclusion of that paper is not changed: as in 0, our gap equation has a nontrivial 
solution at all G > 0, if i7 / 0. Fig. [T] illustrates this important result. Displayed are the r.h.s. 
(curves /) and l.h.s. of equation shown only for the positive valued abscissa, since the curves 
are symmetric. The numerical value of the four-quark coupling is taken to be GA^ = 3. The 
l.h.s. is a constant and indicated by the short-dashed line. The r.h.s. is shown as the long-dashed 
curve (for H = 0) and as full lines for \QH\A~'^ = 0.1 and \QH\A~'^ = 0.5. One sees that in 
absence of the magnetic field one is in the subcritical regime of dynamical symmetry breaking, 
for this choice of GA^. The inclusion of a constant magnetic field, however small it might be, 
changes radically the dynamical symmetry breaking pattern, due to the singular behaviour of 
the r.h.s. close to the origin: the crossing of right and left hand sides of equation ^ does always 
occur and the value of m where this happens is a minimum of the effective potential. In turn 
the trivial solution, which at -ff = is a minimum, is transformed into a maximum for H ^ 0. 

Let us look at these different regimes more closely. As -ff ^ 0, we recover the well-known 
NJL model gap equation 



1 



GN, 
27r2 



\QH\^A^ 



Jo(m2) + 



3m2(A2 -|- m^] 



+ 



(5) 



where Jo{m'^) is 

Jo{m^)=A'^-m^ln(l + ^) . 

\ J 



Equation (l5|) at = admits a nontrivial solution only if r > 1, where r = GA'^Nc/2'k'^ . This 
determines the critical value Gcrit = 2it'^ / h?Nc. The coupling G is said to be over critical if 
G > Gcrit) and subcritical in the opposite case when G < Gcrit- 



^Note that the SU{3) result for the effective potential divided by the number of flavours does not coincide 
with the one flavour case. This is due to group structure factors. 
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0.3 0.4 

mj A 

Figure 1. The the l.h.s. (straight short-dashed hne) and the r.h.s. (other curves) of the gap equation ^ 
related to the four-quark Lagrangian, as functions of mj h.^ all in dimensionless units. Three different 
values of the magnetic field strength H are considered: the long-dashed curve corresponds to if = 0; 
the two full curves represent the cases for magnetic fields of strength jQi/jA^^ = 0.5 (upper curve) and 
\QB.\h.~'^ ~ 0.1. For any H ^ the l.h.s. and r.h.s. intersect and yield the non-trivial solution of ([4]). For 
H = the system is subcritical for the parameter set considered. It becomes critical in the case t > 1, 
where t = GK^N^/iu:^. 



At rr? I <<C. 1 the r.h.s. of equation dH) is 

Here the function t;(^) of the argument = A^/2|(5i?|, is given by 

^;(0 = :3^[l-21nr(e + l)]+V(0- 

This is a monotonically increasing function on the interval < ^ < oo; f (^) = at the point 
^ ~ 1.12; the asymptotic behaviour is 



2? 



1 



2e' 



(7) 
(8) 



where 7 ~ 0.577 is the Euler's constant. 

In the approximation considered the solution of equation (jj]) reads 



Mdyn = Aexp 



To discuss the physical content of this result, we recall that the energy spectrum of relativistic 
fermions in a constant magnetic field H contains Landau levels 



En{Vz) = ±^/rh? + 2\QH\n + pl, 



n = 0,1,2,... 



with pz denoting the projection of the 3-momentum on the z-axis, i.e., along the magnetic field. 
If the fermion mass m goes to zero, as in the present case, there is no energy gap between the 
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vacuum and the LLL. Thus the integer part of + 1 gives approximately the number of Landau 
levels taken into account. 

The first term in equation ([6]) has a clearly defined two-dimensional character with a loga- 
rithmic dependence on the cutoff in the corresponding gap equation 

and, therefore, in the condensate (compare with equation ([5])). Such a behaviour is associated 
with the (1 + l)-dimensional dynamics of the fermion pairing on the energy surface Eq = oi 
the LLL [6]. As long as this term dominates over the second term, v{S,) in ([6]), one concludes 
that the condensate is mainly located on the LLL. Actually this condition is fulfilled nearly 
everywhere at r < 1. This is obvious for = 1, since v{l) = 1/2 -|- = 1/2 — 7 ~ —0.08 
is small compared with l/r. For ^ < 1 we come to the same conclusion after considering the 
asymptotics of the second term ([5]). The other formula, ([7]), can be used to show that the above 
statement is also true for ^ > 1, except near the critical region r ^ 1 — 0, where v{(,) dominates; 
then we must conclude that the condensate spreads over many Landau levels. 

In this special case it is possible to find an analytical solution. Indeed, using d?]) in equa- 
tion ([6]) we obtain 

We further suppose that the two following small variables are of the same order 
m2 f\QH\y 

Then it follows immediately that the term with the logarithm is of order ^/e In y/e and goes to 
zero, when e — > 0. Thus, the gap equation 

is valid only in the region near the critical value r — > 1 — 0. The closer r to 1, the smaller 
is e; Landau levels approach a continuum distribution, and a condensate occupies many levels. 
The physical reason for the changes found in the behaviour of the condensate is the strength 
of the four-fermion interaction which becomes essentially important here. The corresponding 
solution is 

A2 /I 



,^2 \QH\ 
Mdjn = exp 



1 



The near-critical regime which we find here differs from the result of paper [6]. In our case 
the dynamics of the fermion pairing is essentially a (3 -|- l)-dimensional one with a quadratic 
dependence on the cutoff, to be compared with the previous case in equation ([9]), where we had 
a (1 + l)-dimensional behaviour with only a logarithmic dependence on the cutoff. 

3.2 The extended multi-fermion system in a constant magnetic field 

Now we consider the three flavour case with K,gi,g2 7^ 0. In the simplest case with the octet 
flavour symmetry, where current quarks have equal masses rhu = rhd = m^, which we set again 
zero, the system ([2]) reduces to a cubic equation with respect to h = hu = ha = hs 
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0.2 0.4 

m/A 

Figure 2. The effective potential Vh~^ (in dimensionless units) for the parameters GA^ = 3, 
kA^ = -800, AA* = 1.667 • 10^. From bottom to top the curves correspond to \QH\k^'^ = 0.05; 
0.1: 0.2275; 0.3; 0.4. 



with A = gi + (2/3) (72- Making the replacement h = h — k/(36A), one obtains from (llOp 
+ th = b, 



(11) 



where 



4 

t = - 
3 



A V24Ay 



K 


'G 


2 




27A 


X 


-3 


\2A\) 



Am 
3A ■ 



This cubic equation has one real root, if t > 0, i.e., 
A ^ V24A; ' 



(12) 



which is a particular case of the inequalities Assuming that the couplings fulfill condi- 
tion (fT2]) . we find (at any given value of b) the real solution of equation ([TT]) 



h{m) 



D 



with the important property /i(0) = k/36A which ensures that the trivial solution of the gap 
equation, m = 0, exists (see [19] for more details about this part of the potential). 

In Fig. [2] the effective potential is shown for the parameters GA^ = 3, kA^ = —800, AA^ = 
1.67-10^. The catalyzed minima close to the origin at small values of the magnetic field get 
washed out with increasing strength of H giving way to a second minimum at larger values 
of Mdyn- For the lower values of the catalyzed minimum is so close to the origin that one can 
not discern it from the trivial maximum which always occurs at the origin (see also discussion 
around equation (jl4p ). The catalysis of dynamical symmetry breaking can again most easily 
be seen by considering the non trivial solutions of the corresponding gap equation, to which we 
turn now. 

In a constant magnetic field the gap equation is determined by 



27r^h{m) 
h?Ncm 



fim^;A,\QH\). 



(13) 



Comparing this result with equation (jj]), one sees that only the l.h.s. is changed. The six- 
and eight-quark interactions have modified it in such a way that now we get a function h{m)/m 
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u, A # 




Figure 3. The l.h.s. (short dashed line, label u) and r.h.s. (label /) of the gap equation (all in di- 
mensionless units) with the parameter set G, k, A of Fig. [2l From bottom to top the r.h.s. curves have 
|Qi?|A-2 = (long dashed curve), \QH\A-^ = 0.02; 0.1; 0.2275; 0.4 (solid curves). The intersecting 
solutions Afdyn are displayed in Fig. [D 

instead of the former constant term involving only the coupling of four-quark interactions, —1/G. 
The l.h.s. of (113p . abbreviated by u in Fig.O has now a bell-shaped form (short-dashed line), to 
be compared with the horizontal line of Fig. [TJ Since h{m) /m = —1/G + 0{m), the bell-shaped 
curve crosses the ordinate axis at the same point as the former straight line (for the same value 
of GA^). The r.h.s. is again represented by the long-dashed curve {H = case) and by the full 
lines for the finite H cases. As mentioned before these are not altered by the couplings G, k, A. 
The intersection points of the l.h.s. with the r.h.s. curves yield the nontrivial solutions of the 
gap equation: either one or three solutions can be found for m > 0. For the parameter set 
considered the region of coexisting solutions is very narrow, see Fig. [H If equation ()13p has no 
non-trivial solutions at H = (the present case) we say that the set of couplings G, k, A is 
subcritical. It is said to be overcritical in the opposite case. Note that the overcritical set may 
contain G < Ga-it (a parameter set for this case can be found in [22]; there the coexistence of 
solutions fills a larger interval in H, and starts at H = 0). 

The trivial solution, m = 0, corresponds to the point where the potential, V{m), reaches its 
local maximum. The second derivative 



lim 



m) , NJQHl ^ \m\A 

, t; — = lim K — In — — — r 

dm^ m^o 27r2 2\QH\ 



-oo 



(14) 



is negative here. This is the general mathematical reason for the phenomena known as magnetic 
catalysis of dynamical flavour symmetry breaking. The logarithmically divergent negative result 
ensures that this phase transition always takes place, if H ^ 0. This transition does not depend 
on the details related with the multi-quark dynamics, i.e., the result is true even for free fermions 
in a constant magnetic field. 

What is really dependent on the multi-quark dynamics is the local minima structure of the 
theory. Let us recall that in the theory with just four-fermion interactions the effective poten- 
tial has only one minimum at m > 0, and this property does not depend on the strength of the 
field H. We have demonstrated this in Fig. [TJ In the theory with four-, six-, and eight-quark 
interactions one can find either one or two local minima at m > 0. The result depends on the 
strength of the magnetic field H, and couplings G, k, A. This is illustrated in Figs.[2Hll Namely, 
the upper full curve (r.h.s. of equation (fT3j) for \QH\A~'^ = 0.4) has only one intersection point 
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0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 

\QH\/A' 

Figure 4. The non-trivial solutions Mdyn/A of the gap equation, as functions of jQiJjA"^, for the 
parameter set G, k, A of Figs. [2] and [31 For < \QH\A~'^ < 0.2274, Mdyn corresponds to the catalysis 
phase. Then a rapid transition leads the system to the larger Mdyn solution induced by the higher order 
multi-quark forces. A coexistence of both kinds of solutions (at the interrupted line) happens in a very 
narrow window of \QH\A^'^ values (see box insertion: short dashed line are maxima, solid lines minima). 
All quantities are dimensionless. 

with the bell-shaped curve u (l.h.s. of equation ([13]) for GA^ = 3, kA^ = -800, AA^ = 1667.), as 
seen in Fig.[3l This point corresponds to a single vacuum state of the theory at a large Mjyn. As 
opposed to this the curve with \QH\A~'^ = 0.2275 has three intersections with the same curve u; 
these intersections, taken in order, correspond to a local minimum, a local maximum and a fur- 
ther local minimum of the potential (see also Fig. [4]) . Reducing further the field strength one ob- 
tains again only one solution, at considerable smaller values of Mdyn, corresponding to catalysis. 

Therefore the first minimum catalyzed by a constant magnetic field is then smoothed out 
with increasing H. It ceases to exist at some critical value of \QH\A~'^, from which on only 
the large Mdyn solution survives, as shown in Fig. HI for the parameter set of Figs. [21 [3l This 
process is accompanied by a sharp increase in depth of the effective potential at the second 
minimum (see Fig. [2]). It is clear that the only way to wash out the first minimum (due to 
equation ()14p ). is by lowering the barrier between this state and the second minimum. We 
also observe that the second minimum is unremovable, because the asymptotic behaviour of 
the functions in equation ()13p is such that the l.h.s. dominates over the r.h.s. at large m/A. 
Let us stress that this scenario will always be possible if at = the six- and eight-quark 
interactions induce dynamical symmetry breaking in the subcritical regime G < Gcrit of four- 
quark interactions (the case presented in [22]). If the system remains subcritical also for the 
higher-order interactions, as in the example shown here, the regime of coexistence of more than 
one minimum of the gap equation at finite H shrinks eventually to zero by lowering further the 
eight-quark interaction strength A (keeping the remaining parameters fixed). Then one can still 
observe a crossover transition from small to large Mdyn solutions. For example by changing the 
r.h.s. of the stability condition G — (^)^ I ^ ^ from 2.33 (parameter set shown in Figs. [2H1]) 
to 1., one obtains qualitatively the same picture as in Fig. [H (of course without coexistence 
region) but a smoother transition. 

4 Concluding remarks 

The effects of interactions among fermions on catalysis of dynamical chiral symmetry breaking 
by a constant magnetic field have been considered. A generalized NJL model with U{3) flavor 
symmetry containing up to eight-quark interactions has been at the basis of this study, and the 
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motivation for the use of such a Lagrangian was briefly reviewed. A critical survey of the non- 
trivial solutions of the gap equations has been obtained and compared to the known solutions 
of the conventional four-quark interaction NJL model. 

We have discussed two cases: (i) all coupling strengths are subcritical, (ii) the addition of six 
and eight quark interactions transform the original subcritical G < Gcrit four-quark system into 
an overcritical one, inducing the formation of a large condensate. 

In both cases we have that in the region m^/A^ <C 1 the four-fermion interaction dominates 
the behaviour of the system. Since their coupling strength is small, G < Gcrit; the massless 
fermions behave like almost free particles moving in a weak external magnetic field, with access 
to a large number of Landau levels, ^ ^ 1. This field catalyzes the process of fermion-antifermion 
pairing on the energy surface E'o = of the LLL. The first minimum localized at rr? j ^ 1 is 
exactly formed by such an (1 -|- l)-dimensional condensate. If six- and eight-fermion forces would 
not act on the system, this ground state would be stable: our formulae as well as the result of 
paper [6] show clearly that a slow increase of the strength H does not wash out the condensate 
from the energy surface E'o = of the LLL. However, when the six- and eight-fermion interactions 
are present, a slowly increasing magnetic field destroys finally this ground state. In case (ii) 
this always happens; in case (i) the transition can either be sharp as in (ii) or of a crossover 
type, depending on the interaction strengths. The new condensate has a (3 -|- l)-dimensional 
structure similar in every respect to the standard NJL case with broken chiral symmetry at 
H = 0, i.e., when the condensate spreads over many single fermion states. This is because the 
increasing magnetic field enlarges the dynamical fermion mass, and scales of order m/A ~ 1 
become relevant. At these scales the 't Hooft and eight-quark interactions push the system to a 
new regime, where the fermions are not anymore free-like particles: they interact strongly with 
each other and this interaction changes essentially the fermionic spectrum and the structure of 
the ground state with all the above mentioned consequences. 

Thus we have obtained not only a correct description of the well-known physics related with 
the LLL, but have found also a clear signature for the possibly important role played by 't Hooft 
and eight-quark interactions. Namely, in the presence of these interactions the magnetic field 
can change the condensation zone from the zero-energy surface of the LLL to a wide region 
spread over many Landau levels and vice versa. One can expect that hard gamma emissions 
accompany this process. 

The value of the characteristic scale A has not been specified yet. We assume that this value 
is determined by the problem under consideration. Its choice can also be motivated by the 
number of Landau levels to be considered. The characteristic scale of the magnetic fields which 
can induce such a transition is of the order H = 7 • 10^^ A^ Gauss/MeV^s, and actually depend 
on the cutoff involved in the problem. For instance, in hadronic matter it is probably reasonable 
to assume that A ~ 800 MeV, leading to H = 4.5 • 10^^ G. 

A potentially interesting area where this effect may find applications are studies of com- 
pact stellar objects in presence of strong magnetic fields, in particular the young neutron stars, 
magnetars [M]- The surface magnetic fields are observed to be > 10^^ G, but actually they 
can be even much higher at the core region. Another area is connected with the electroweak 
phase transition in the early Universe [35], where the strength of magnetic fields can reach 
H ~ 10^"^ G. Even giving conservative limitations on the primordial magnetic field of present 
epoch value of 10~^ G, based on primordial nuclear synthesis calculations [36], it seems clear 
that the magnetic field reached the transition regime discussed in this paper in the past. A very 
rough estimate based on the relation H ccT^ {Ty being the neutrino temperature), gives a mag- 
netic field strength at baryogenesis of at least 10^^ G, of 10^^ G at the epoch of electroweak 
breaking and of 10^^ G at the GUT scale. The coupling to those strong magnetic fields will 
have influence on those processes. However to incorporate the findings of the present work into 
those mechanisms is beyond the scope of the present article. 
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